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139 La nuclear quadrupole resonance measurements in lightly doped L^Cui^LLCU have been per- 
formed to reveal the dependence of the magnetic properties of the antiferromagnetic CuCh planes 
on the character of the doped holes and their interactions with the dopant. A detailed study shows 
that the magnetic properties are remarkably insensitive to the character of the dopant impurity. 
This indicates that the added holes form previously unrecognized collective structures. 
PACS numbers: 75.30.Kz, 76.60. Jx, 74.72.Dn, 76.60. Gv 



Full understanding of the character of holes added 
to cuprate planes and their interactions with the two- 
dimensional lattice of Cu spins remains a crucial and un- 
solved problem in the high temperature superconductors. 
While the detailed mechanism is poorly understood, the 
rapid suppression of the antiferromagnetic (AF) ordering 
temperature Tv by doping is clearly related to the dis- 
ruptive effects of mobile holes: < 3% added holes whether 
from Sr substitution, addition of interstitial oxygen or in- 
plane substitution of Li for Cu |Q suppresses Tv to zero, 
yet ~ 30% isovalent substitution of Zn or Mg for Cu is 
required to produce the same effect. Li and Sr-doped 
holes have very different mobilities. For x or y ss 0.025, 
the room temperature resistivity p of La2Cu 1 _ 2 .Li 2; 04 
exceeds that of La2-j / Sr y Cu04 by over an order of mag- 
nitude mm ; more strikingly, for Sr-doping dp/dT > 
for T > 100 K, in contrast to the negative slope found in 
Li-doped material for all x and T. 

It is well recognized that the 2D cuprates are inclined 
toward microscopic charge inhomogeneity [^-Q. Evi- 
dence for such an effect in lightly doped La2-j,Sr a CuC>4 
was obtained from a scaling analysis of the doping y and 
temperature T dependence of the static susceptibility 
x{y,T) [pj which indicated that the magnetic correlation 
length is limited to the dimensions of AF domains (finite- 
size scaling) formed by microsegregation of doped holes 
into hole-rich domain walls surrounding hole-free, AF do- 
mains. Interpretations involving charge-stripes have also 
been proposed ||. Castro Neto and Hone |l(J have ex- 
amined the influence of doping on the long wavelength 
properties of a 2D antiferromagnet in a model in which 
charged stripes cause the exchange coupling J to become 
anisotropic; this model reproduces the relationship be- 
tween M° (M s is the sublattice magnetization; M s ° is 
that obtained by extrapolation of data for T > 30 K to 
T = 0) and Tv as the two are suppressed by Sr-doping 
in La2_j,Srj,Cu04 JlT| . However, using a similar model, 
van Duin and Zaancn (l2) find that Tjv is suppressed 
much more rapidly than M s ° with increasing anisotropy 
(doping). 

We have used 139 La nuclear quadrupole resonance 



(NQR) measurements to microscopically examine the 
effects of doped holes on the AF spin correlations in 
La 2 Cu 1 _-_i_0 4 (0.019 < x < 0.025). We find that the 
magnetic behavior of lanthanum cuprate is remarkably 
insensitive to the detailed nature of the dopant, in spite 
of the differing charge transport associated with the two 
dopants. In addition to the similarly strong suppression 
of Tjv by doping, we find the identical correspondence 
between the suppression of M° and Tv by doping which 
has been observed in La 2 -, y Sr, / Cu04. Further we show, 
for the first time, that in the vicinity of Tv the dynamical 
susceptibility, as reflected in the nuclear spin relaxation 
rate 2W, follows a scaling law consistent with the finite- 
size scaling demonstrated in the static susceptibility by 
Cho et al. [B. Finally, at low temperature, we find that 
two peculiar features are very similar in the two systems. 
These are the very strong peak in 2W (at a temperature 
Tf = 10-16 K depending on x) that indicates freezing of 
spin-degrees of freedom, which is accompanied at slightly 
higher temperatures (~ 30 K) by the abrupt recovery of 
M S (T), almost to x = values. 

Thus, while very small concentrations of added holes 
induce a range of characteristic magnetic properties 
which are entirely insensitive to the nature of the dopant, 
the transport properties are very sensitive to the dopant. 
Unable to understand these contrasting behaviors as aris- 
ing from properties of individual holes, we conclude holes 
form collective structures. We will argue that holes 
form charged, anti-phase domain walls |13j which sur- 
round mobile domains in which the phase of the AF or- 
der is reversed. Such mobile domains will suppress the 
time-averaged static moment thus suppressing AF order 
and M s . These domain structures will have contrast- 
ing interactions with in-plane vs. out-of-plane dopants 
(e.g., stronger scattering by in-plane impurities) which 
explain the different transport behaviors, while the uni- 
versal magnetic properties can be understood as long as 
the domains are sufficiently mobile that they move across 
a given site rapidly compared to a measurement time. 

Three powder samples of La2Cu 1 _ 2 .Li x 04 (labeled Al, 
Bl and B2) were prepared from starting material contain- 



ing concentrations x nom of Li as described elsewhere . 
Table || shows various measured properties of these sam- 
ples; definitions of the parameters shown and the means 
by which they were determined will be discussed in what 
follows. All the samples were annealed in flowing nitro- 
gen gas after sintering to minimize excess oxygen, whose 
motion can contribute to 139 La spin-lattice relaxation at 
high T @Q. Annealing in nitrogen eliminated this con- 
tribution in our samples. The smaller value of x relative 
to x nom for sample Al was deduced from its higher value 
of T N relative to sample Bl. 139 La (/ = |) NQR and re- 
laxation rate measurements in La2Cu 1 _ 2 ,Li a; 04 were per- 
formed for 4 K < T < 300 K on both the 2v Q (±| <-»±§) 
and 3vq (±|<->±|) transitions. The spectra were ob- 
tained by plotting the integrated intensity of the spin- 
echo signal as a function of spectrometer frequency. 
The nuclear spin-lattice relaxation rate was measured by 
monitoring the recovery of the magnetization after satu- 
ration with a single ^ pulse. 

The 139 La NQR results for the three samples are sum- 
marized in Fig. 0. Below Tjv, the ordered Cu moment 
generates an internal magnetic field H at the 139 La site 
which splits the NQR line into a doublet with frequencies 
v\ and V2- The magnitude of the splitting, A = v\ — V2, 
is a direct measure of the component H z — H cos 9 of the 
internal field H at the La site along the principal (z) axis 
of the electric field gradient (EFG) tensor |9j (9 is the an- 
gle between H and z). The two frequencies for the 2vq 
transition are given by v\ 2 = 2vq ± {^ n /2'K)H z where 
the nuclear gyromagnetic ratio 7 Tl /27r = 601.44 Hz/G. 
As shown in Fig. |](a), for T > 30 K, A is suppressed 
by doping, and then displays an abrupt increase below 
30 K. Similar effects have been observed in recent /iSR 
and neutron scattering measurements in La2Cu 1 _ a ,Li :I ;04 
fl5f . The solid curves in Fig. 0(a) are fits of the temper- 
ature dependence A(T) = A (l - T/Tn) 13 to the data 
for T > 30 K; we find f3 = 0.44 ± 0.01. The magnitude 
and temperature dependence of vq for a given doping is 
almost identical to that found in Ref . [[fl| ; the linewidths 
are also quite comparable: 110, 120, and 140 kHz at 250 
K for our samples Al, Bl and B2, respectively. 

Fig. |l|(b) shows the T-dependence of 2W , which dis- 
plays a strong peak at T — Tf — 10-16 K and a weak 
peak in the vicinity of Tat. 2W was obtained by fitting 
the recovery data to the theoretical expression for mag- 
netic relaxation and for single pulse saturation For 
T > Tf, this expression fits the data well, indicating a 
single rate arising from magnetic fluctuations. Although 
the rate becomes distributed below Tf (consistent with 
freezing in an inhomogeneous distribution of internal 

TABLE I. Properties of the three samples: Al, Bl and B2. 
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0.020 


0.019 
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0.042 


1.52 


Bl 


0.020 


0.020 


145±5 


123 


193 


15 


163 


0.038 


1.52 


B2 


0.025 


0.025 


85±10 


117 


120 


16 


114 


0.020 


1.25 



fields), the same fitting procedure was applied to data 
for the first decade of recovery. While this increases the 
uncertainty in 2W, we find that varying the fitting proce- 
dure has essentially no effect on the position of the peak 
at Tf. 

As seen in Fig. 0(b), 2W(T) is the same at the 2vq 
and ?>vq transitions in sample Al, establishing that 2W 
is due to a magnetic rather than structural mechanism. 
In addition, we note that 2W decreases with increasing T 
for T > Tn in contrast with the results in La2_ l/ Sr !y Cu04 
which show an enhancement of 2W at high T due to 
the motion of excess oxygen Jl4|| . Thus, annealing to 
remove excess oxygen enables us to obtain, for the first 
time, intrinsic data for spin dynamics above XV that are 
isolated from any significant influence by mobile oxygen. 

Our discussion will focus on two aspects of the data: 
(i) the anomalous behavior of both the static and the 
dynamical magnetic properties at low T, and (ii) the 
broad and weak peak in 2W around Tjy, which contrasts 
with the sharp peak at Tjv observed in, e.g., undoped 
Sr2Cu02Cl2 These results are, overall, very similar 
to those found in La2- y Sr y Cu04 0JllJ|. 

The strong low-T peak in 2W clearly indicates freezing 
of spin degrees of freedom. Analyzing the data in terms 
of activated behavior, [2H^(T) oc exp(EafkBT)} as shown 
in Fig. ||, gives values of E a (see Table |) similar to those 
in La2-j,Srj ) Cu04; where Ea/ksTf = 8.9, and Tf ~ 11 
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FIG. 1. " y La NQR in LaaCui^Li^CU: 
(a) A = v\ — V2 vs. T. Solid curves are fits to the behav- 
ior A(T) = A (l -T/T N ) 9 , (3 = 0.44±0.01. The inset shows 
a representative NQR spectrum at 4 K. (b) 2W vs. T. 
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-16 K E a ~ 100-145 K 0. Thus, the spin freez- 
ing is independent of the dopant, and, consequently, spin 
freezing is not solely determined by binding of holes to 
the dopant. We note that Tf{x) (Table does not sat- 
isfy the empirical relation Tf = (815 K)x obtained in 
Rcf . [[llj , but the relationship between Tf and T/v is the 
same in both systems. Such disagreement may arise from 
a discrepancy in actual hole-doping level. 

We turn now to the behavior of 2W in the vicinity of 
Tjy. We find a strong doping dependence of the width of 
the 2W peak as found in the static susceptibility by Cho 
et al. j^], and our results support the finite-size scaling 
proposed there. To isolate the behavior near Tjv from 
contamination by the tail of the spin-freezing peak at 
Tf, we first correct the 2W data by subtracting the fit- 
ting results in Fig. || (solid lines). Results are shown in 
Fig. ||(a). The observed behavior cannot be understood 
in terms of conventional critical behavior |l7|| . There, the 
T-dependent growth of £ is set by the exchange coupling 
constant J, and this determines the width of the peak on 
the high-T side. The strong x-dependence of the width 
then leads to the implausible conclusion that J is strongly 
x-dependent. Nonetheless, it is clear that the behavior 
of 2W around Tjv should reflect cooperative behavior of 
correlated spins near their ordering temperature. 

The peak in the static susceptibility observed at Tm is 
due to the small canting of the ordered moment which 
arises through the Dzyaloshinski-Moriya interaction. In 
undoped La2Cu04 the shape of this peak is well de- 
scribed by theory fL8|| , however with doping the peak 
broadens and this description is no longer adequate. 
Cho et al. argued that rather than increasing exponen- 
tially with decreasing temperature, £ (T) in lightly doped 
La2-j,Sr a Cu04 is limited by confinement to AF domains 
of size L defined by hole- rich domain walls Q . They were 
successful in describing the doping dependent shape of 
the peak in x(y> T) (for T > Tjv) for a range of y through 
the scaling relation x{v, T ) = x{f{y)[ T - T N{y))}, and 
they identify f(y) ~ 0.02/y 2 with L 2 . If this peak shape 
is indeed due to £(T), consistency of the scaling of 2W/T 
and the static susceptibility is expected since both the 
static and dynamic susceptibilities (\' and \" > respec- 
tively) are determined by the same T-dependent £. 
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FIG. 2. 2W vs. 10/T. Solid lines are fits to the activated 
behavior 2W oc exp(E a /ltBT). 



different from the peak in x(T) [marked by arrows in 
Fig. ||(a); see also Table [|. Given the incomplete under- 
standing of x(T,x) for finite x, the peak in 2W is likely 
the better indicator of the temperature at which local AF 
ordering occurs. Although not understood, we believe 
this phenomenon is a general feature of lightly doped 
lanthanum cuprate, and is more prominent in our data 
because the obscuring effects of mobile excess oxygen at 
temperatures comparable to Tjv have been eliminated 
here. We find that 2W(T, x) exhibits the same scaling be- 
havior with the same scaling function f(x) = (0.02/x) 2 . 
The strength of the peak in 2W is strongly x-dependent, 
and we apply a normalization factor C(x): 



^{x,T} = C(x)^{f(x)[T- 



n(x)}}. 



(i) 



The data scaled in this way are shown in Fig. [|(b), and 
we find C{x) ~ (.004/x) 2 as shown in Table [j| Given the 
discrepancy between Tpj and T^, these data cannot be 
taken as demonstrating scaling of 2W; however, our find- 
ing that the same scaling function fix), when combined 
with a normalization factor also proportional to 1/x 2 , 
leads to consistent scaling behavior provides important 
corroboration for the finite-size scaling hypothesis of Cho 
et al §. 

The reduction of 2W/T associated with limitation 
of the size of correlated domains is not unexpected. 
The spin-lattice relaxation rate is given by 2W/T ~ 
J2 q jS ^{q)[x"{qi UJ )/m\, where A(q) is the (sample inde- 
pendent) hyperfine coupling constant. If the spectrum of 
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FIG. 3. (a) T— dependence of 2W: 2W data are corrected 
for the contribution from the spin-freezing effects represented 
by the solid lines in Fig 2. The arrows indicate Tjv- (b) Plots 
of normalized 2W/T as a function of f(x)(T — T^) with 
f(x) = (0.02/a;) 2 (see text for details). 
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spin fluctuations is described by a Lorentzian distribu- 
tion with a characteristic frequency T, then x"{Qi ~ 
x'(3)/r; this highlights the additional sensitivity of 2W 
to the characteristic fluctuation rate Y compared to the 
real part of the susceptibility. As the number of spins in a 
domain (argued in the finite-size hypothesis to go as 1 jx 2 
[H) increases, Y will decrease; and, thus, the relaxation 
rate will increase with decreasing doping. Qualitatively, 
this is observed: 2W/T ~ C(x) ~ (l/x) 2 . 

The detailed similarity of the magnetic properties of 
Sr and Li-doped lanthanum cuprate cannot be under- 
stood in terms of individual holes interacting with impu- 
rities; this suggests that holes form collective structures. 
A full explanation of these magnetic properties cannot 
be given at this point, but to be concrete we propose a 
picture that can successfully account for the distinctive 
and complex behavior of A and 2W below Tjv. Phase 
segregation of holes to boundaries of hole free regions ||] 
or to hole rich "rivers" Jj| has been proposed. Based 
on the results of Tranquada et al. [ p"3| at higher dop- 
ing in Lai.48Ndo.4Sro.i2CuC>4, we suggest a contrasting 
picture in which the stripes are densely populated with 
holes (e.g., 1/2 hole per stripe Cu site) and constitute 
anti-phase domain walls. We further suggest that these 
stripes close to form mobile loops that surround small 
domains (which we call anti-phase domains) in the anti- 
ferromagnet. Such loops have been theoretically argued 
to be the stable configuration at low doping jl9| . Because 
these encircling stripes are anti-phase domain walls, the 
phase of the AF correlations in an antiphase domain is 
reversed with respect to the background. Motion of these 
domains (above 30 K) over a particular site thus reverses 
the spin orientation and reduces the time-averaged static 
moment at that site. 

To estimate the size and density of domains required 
to explain experimental data, we consider a simple model 
pof in which there is, on average, one domain for ev- 
ery N Cu sites (separated by lattice constant a). These 
N sites occupy an area L 2 = Na 2 , and, of these, 
N- sites comprise the anti-phase domain whose width 
is I = a^J N-. The fractional reduction of the time- 
averaged static moment at a given site will then be 
(N + - N-)/N where N = N + + N-. We set this 
equal to R(x) = M°(x)/M°(x = 0). Assuming 1/2 hole 
per Cu in the stripe, we have 4Z ~ 2xNa (we assume 
that all donated holes are incorporated into loops), and 
I /a ~ (1 — R)/x. Using measured values of R(y) for 
La2- J/ Srj / Cu04 |)|, we find l/a increases with doping to 
~ 20 at y = 0.018 where R has decreased to ~ 0.6. Re- 
covery of the sublattice magnetization below 30 K occurs 
when the domains either become pinned to the lattice 
or evaporate as the constituent holes become pinned to 
donor impurities. Mobile anti-phase domains will intro- 
duce disorder and reduce the total inter-plane coupling 
of AF ordered regions thus reducing Tv Q ■ 

In conclusion, we have presented a complete set of 
139 La NQR data in lanthanum cuprate doped by in-plane 
substitution of Li for Cu. The most striking result is 



the remarkable insensitivity of the magnetic properties 
of the AF Cu02 planes to the nature and location of the 
dopant; this in spite of marked differences in transport 
properties. This suggests the magnetic properties reflect 
collective hole phenomena. The behavior of 2W around 
Tjv cannot be understood in the context of conventional 
critical behavior; we have demonstrated consistency with 
the finite-size scaling approach proposed earlier ||. We 
propose that all of the data can be understood by means 
of a model in which doped holes segregate to domain 
walls that enclose mobile, anti-phase bubbles which re- 
duce the time-averaged ordered moment thus suppressing 
M° and Tv. 
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